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Abstract 

We perform a molecular dynamical study of the isolated d = 1 classical 
Hamiltonian H = \ J2iLi Lf + J2i=£j 1 ~ cos ^ t ~ 9 ^ ; (a > 0), known to exhibit a 
second order phase transition, being disordered for u = U/NN > u c (a, d) and 
ordered otherwise (U = total energy and = — \l a J^ d )- We focus on the 
nonextensive case a/d < 1 and observe that, for u < u c , a basin of attraction 
exists for the initial conditions for which the system quickly relaxes onto a 
longstanding metastable state (whose duration presumably diverges with N 
like N) which eventually crosses over to the microcanonical Boltzmann-Gibbs 
stable state. The temperature associated with the (scaled) average kinetic 
energy per particle is lower in the metastable state than in the stable one. It 
is exhibited for the first time that the appropriately scaled maximal Lyapunov 
exponent (metastable) oc N~ Kmetastable ; (N — * oo), where, for all values 
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of a/d, K me tastabie numerically coincides with one third of its value for u > u c , 
hence decreases from 1/9 to zero when a/d increases from zero to unity, re- 
maining zero thereafter. This new and simple connection between anomalies 
above and below the critical point reinforces the nonextensive universality sce- 
nario. 

PACS numbers: 05.70.Fh, 64.60.Fr, 05.50.+q 

The foundations of statistical mechanics, hence of thermodynamics, is a subtle and fas- 
cinating matter which has driven enriching controversies and clarifications since more than 
one century (see, for instance, Einstein's remark on the Boltzmann principle [0]). The field 
remains open to new aspects and proposals. One of these is nonextensive statistical mechan- 
ics, proposed in 1988 (see || for reviews). This formalism is based on an entropic index 
q (which recovers usual statistical mechanics for q = 1), and has been applied to a variety of 
systems, covering certain classes of both (meta)equilibrium and nonequilibrium phenomena, 
e.g., turbulence [|J, hadronic jets produced by electron-positron annihilation ||, cosmic rays 
H, motion of Hydra viridissima J7J, among others. In addition to this, it has been advanced 
that it could be appropriate for handling some aspects of long-range interacting Hamiltonian 
systems. This possibility is gaining plausibility nowadays, as argued in |§ and elsewhere. 
Indeed, in molecular dynamical approaches of isolated systems, strongly nonmaxwellian ve- 
locity distributions have recently been observed that are consistent with such possibility || . 
A paradigmatic system in the realm of this discussion is the following classical Hamiltonian: 

Z i=l Vftj 'ij 

The inertial planar rotators (ferromagnetic XF-like model) are localized at the sites of a d- 
dimensional periodic lattice. As distance (measured in crystal units) for a given pair (i, j) 
we consider the shortest among all the possible ones (due to periodicity). For d = 1 we have 
rij = 1, 2, 3, for d = 2 we have r^- = 1, a/2, 2, for d = 3 we have = 1, v2, \/3, 2, 
and so on. The so called HMF system is recovered for a/d = 0, and the first-neighbor 
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model is recovered for a/d — > oo. Hamiltonian (1) is extensive if a/d > 1 and nonextensive 
if < a/d < 1. This can be seen as follows. If we define 

r Nl/d , , N~L-a/d — q, I A 

N=l+d drr d - 1 r- a = - (2) 

Ji 1 — a/d 

it can be easily checked that the energy scales as NN, i.e., it is asymptotically proportional 
to N if a/d > 1, to TVlniV if a/d = 1, and to N 2 ~ a / d if < a/d < 1. This Hamiltonian is 
sometimes presented in the literature in the following form: 

n , = 1 ^f + W 1 - c °f- 9i) <«*<>). (3) 

Z i=l iV i# 'ij 

which artificially makes its energy to scale as N, W(a/d). The transformation from this form 
to the one presented in Eq. (1), adopted from now on in the present work, has been described 
in detail in |1(J. This system has since long been shown JTT ] to obey Boltzmann-Gibbs (BG) 



statistical mechanics for a/d > 1. What happens for < a/d < 1 is a subtle question 
which is under intensive study nowadays |T2|-|iT||. In fact, several long-range-interacting 



systems are since long known fll8|-c22| to present a variety of thermodynamical anomalies, 
such as negative specific heat and super diffusion among others. The molecular dynamics 
in the isolated Hamiltonian (1) with total energy U exhibits, for infinitely large time, the 
existence of a second order phase transition at u = U /NN = u c (a, d). For u>u c the system 
is disordered (paramagnetic-like); otherwise, it is ordered (ferromagnetic-like). It exhibits 
anomalies on both sides of the critical point. 

For u > u c (i.e., in the disordered phase), after a quick transient, the one-particle dis- 
tribution of velocities gradually becomes Maxwellian in the iV — >• oo limit, in accordance to 
what is expected within BG statistical mechanics. However, while N increases, the entire 
Lyapunov spectrum approaches zero, which is a quite anomalous behavior; indeed, no such 
weakening of chaos is expected nor observed for a/d > 1. This weakening of the sensitivity 
can be characterized through the maximal Lyapunov exponent A™^ c (appropriately scaled 
as indicated in JTUJ) which, in the N — ► oo limit, vanishes as oc N~ Kd (d stands for 

disordered phase); decreases froml/3 to zero while a/d increases from zero to unity, and 



remains zero thereafter WM. The fact that the sensitivity to the initial conditions becomes 



sub-exponential (possibly a power-law) strongly reminds what has been observed [p3|-|3T 
in a variety of low dimensional maps, which are known to be adequately described within 
nonextensive statistical mechanical concepts. 

For u < u c , after a quick transient, the behavior depends from the initial conditions. 
Two wide basins of attraction exist in the space of the initial conditions. One of them 
(which includes Maxwellian velocity distribution and all angles equal) yields a standard BG 
microcanonical distribution which approaches the BG canonical one in the limit iV — > oo. 
The other one (which includes waterbag and double waterbag velocity distribution and 
all angles equal) yields a longstanding metastable (quasi-stationary) state (whose associated 
magnetization is basically zero) and only at very large time joins the BG distribution (whose 
associated magnetization is nonzero). The duration r of this metaequilibrium state diverges 
with N. It has been conjectured [^2| that it does so as r oc N. Recent results support this 
scaling; indeed, (i) for a = 0, this conjecture implies r oc N, which has been verified ||, 
(ii) for fixed N, it implies that r exponentially decays with a, which once again has been 
verified KIT 



Our focus in this paper is on the metastable state of the d = 1 model, which we study for 
typical values of (a,u,N). The time evolution of the model has been generated integrating 
the equations of motion through a 4 th order symplectic algorithm [33| with a relative error 
in the total energy conservation less than 10~ 4 . We verify that the time evolution of the 
(scaled) average kinetic energy per particle (which plays the role of temperature) exhibits 
two plateaux, the first one being anomalous and the second one being of the BG class. This 
is illustrated in Fig. 1 for a = 0.6. In the same figure we show the time evolution of A™<^ c . 
As in the a = case, one expects also for < a/d < 1 two plateaux in AJ^ c (i). We 
can see, however, that for a = 0.6 the difference is almost unperceptively small; it might 
happen that this difference quickly decreases with a, as it is the case for r, but such study 
is out of the scope of the present work. The systematic detection of both plateaux in the 
temperature enabled the calculation of the caloric curves, as illustrated in Fig. 2. We clearly 



see the existence of negative specific heat for the metastable state, just below u c . Then by 
focusing on small time (after the transient nevertheless), it was possible to calculate the 
iV-dependence of {metastable) for typical values of a. The results are shown in Fig. 3. 
We verify that \™<* c (metastable) oc N~ K,meta3table , where n m etastabie decreases from 1/9 (thus 
confirming |34}| ) to zero, while a increases from zero to unity, and remains zero thereafter. 
Furthermore, we numerically verify a remarkable property, namely (see Fig. 4) 

«— «. = y (Va). (4) 

This constitutes the first connection found for this type of models between the anomalies 
below and above the critical point. This is a conceptually important point. Indeed, if 
nonextensive statistical mechanics is relevant for such long-range interacting systems as the 
velocity distributions presented in |8j seem to suggest, one would expect the model to be 
somehow associated with a single value of the entropic index q for all energies, both below 
and above possible critical points. Eq. (4) makes this possibility plausible. Before ending 
let us mention that no anomalies were detected nor expected for X^^stable) (i.e., in the 
BG regime emerging at large time), which should gradually become positive N- independent 
values for all values of a. This is of course consistent with the picture that lim^oo liHijv-^oo 
(anomalous thermodynamical metaequilibrium) and liniAr^oo lim^oo (BG thermodynamical 
equilibrium) are not! interchangeable if < a/d < 1, whereas they are if a/d > 1. 
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CNPq, PRONEX, FAPERJ (Brazilian agencies), and FCT (Portugal) is also acknowledged. 
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Fig. 1 - Time evolution for twice the (scaled) average kinetic energy per particle 



{Ekin)/NN (which plays the role of temperature; upper curve) and the (scaled) largest 
Lyapunov exponent A™<^ c (lower curve) for a = 0.6, u = 1 and iV = 1000. We have aver- 
aged 10 different water-bag initial conditions for the velocities (all angles were initially set 
parallel to each other). 

Fig. 2 - Microcanonical caloric curves for typical values of a and N. The lower branch 
corresponds to the metastable state. The stable state is indicated with the dashed line. 

Fig. 3 - iV-dependance of the largest Lyapunov exponent \™* c (metastable) for a 
ranging from to 1.2. The average in the interval 10 < t < 3000 has been considered as 
the metastable state value (the very slight increase of the Lyapunov exponent occasionally 
observed up to t = 3000 is numerically without consequences). 

Fig. 4 - a/(i-dependance of 3 x K met astabie (full cercles). Open triangles, cercles and 
squares respectively correspond to K d of the d = 1, 2, 3 models fT0| , |14fl . The arrow points to 
1/3, value analytically expected |15 , 16] to be exact for a = and u > u c . 
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